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It is shown that it is possible to perform a model-independent extraction of the CKM Unitarity 
Triangle angle 7 using only the decays i?" — > DK'^-k~ with D K%-k^-k~ and DK'^-k~ 
with flavour-specific D decays. The proposed method can also utilise the B" — > DK^n^ data with 
CP-eigenstate decays of the D meson in a model-independent way. 
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I. INTRODUCTION 

Among the fundamental parameters of the Stan- 
dard Model of particle physics, the angle 7 = 
aig {-VudV*i,/VcdV*ij) of the Unitarity Triangle formed 
from elements of the Cabibbo-Kobayashi-Maskawa quark 
mixing matrix [l|, Q has a particular importance. Not 
only is it one of the least well determined fundamen- 
tal parameters of the Standard Model, it is also the only 
CP-violating parameter that can be measured using only 
tree- level decays 0. The precise determination of 7 is 
thus a critical element of the programme towards under- 
standing the baryon asymmetry of the Universe, and is 
one of the main objectives of planned future B physics 
experiments (see, for example, 04^)- 

A method to measure the CKM phase 7 in 
DK~^TT~ decays by comparing the Dalitz plot distribu- 
tion obtained when the neutral D meson is reconstructed 
in CP eigenstates to that when flavour-specific states are 
used has been proposed recently Q. The method builds 
on the original proposal of Gronau, London and Wyler 
(GLW) for a similar analysis using i3+ DK+ de- 
cays, and exploits the fact that CP violation effects are 
expected to be enhanced in — )■ DK*^ decays since the 
interfering amplitudes are of comparable magnitude [l^l ■ 
Since the charge of the kaon in — >■ DK~^Tr~ decays un- 
ambiguously tags the flavour of the decayin g B meson, 
there is no need for time-dependent analysis 

Com par ed to previous proposals for similar analy- 
ses |13 - [l4| this method has several attractive features: 
1) the amplitudes can be measured directly in the 
Dalitz plot analysis, relative to the flavour-specific B'^ 

amplitude; 2) there is no need to normalise to 
decays with another final state; 3) there is sufficient in- 
formation to extract 7 using only D decays to CP-even 
states; -'^ 4) the sensitivity to 7 docs not depend strongly 
on the values of 7 and strong phases involved; 5) 7 is 
determined with only a single unresolvable ambiguity 



CP-odd D decays with large branching fractions such as D — > 
Ksn^ are difficult to reconstruct, particularly in experiments at 
hadron colliders. 



(7 — > 7 + tt) . These advantages have much in common 
with those of the P+ DK+, D K^tt+tt^ method, 
which is currently providing the most constraining mea- 
surements of 7 Moreover, this technique is espe- 
cially promising for the LHCb experiment since it does 
not require the reconstruction of neutral particles in the 
final state. 

There is, however, a potential drawback of the method, 
in that the fit to the Dalitz plot distribution is inher- 
ently model-dependent. The magnitude of the model- 
dependence has been estimated to be small [13, but in 
view of the aim to achieve an 0(1°) measurement of 7 
at future experiments, it is certainly desirable to explore 
possibilities to remove such systematic uncertainties. 

A model-independent analysis procedure has been pro- 
posed for the extraction of 7 from P+ — DK~^, D — 
A'gTT+TT^ decays [Tsl. [Tgj. The method uses data from 
charm factory experiments - specifically, CP-tagged 
charm mesons - to obtain information about the strong 
phase variation across the D decay Dalitz plot. A sim- 
ilar approach to probe the strong phases present in the 



P" 



P/v+TT Dalitz plot is not immediately available 



since there is no realistic way to obtain sufficient samples 
of CP-tagged B decays. 

In this paper we show that the study of Dalitz plot 
distributions in the decay chain P° DK^tt^, D ^• 
KgTr'^TT~ allows to translate the model-independent mea- 
surement of the D — >■ Kg7T~^7T~ amplitude to the P° — >■ 
DK~^TT~ mode. This approach allows not only to mea- 
sure 7 model-independently in this decay chain, but also 
to use the obtained constraints on the P° — >■ DK^tt" 
amplitude in the P" — > DcpK~^Tr~ Dalitz plot analysis. 
For clarity, our method is independent of model assump- 



tions in both P 







DK^ 



and D 



decays. 



The more straightforward cases which arc either model- 
dependent in P'^ — > DK^TT^ but model-independent in 
D — > /\^7r+7r^, or vice- versa, are also possible but are 
not discussed in detail. A fully model-dependent analysis 
is, of course, also possible. 

The remainder of the paper is organised as follows. In 
Section [IT] we outline the basic idea and the formalism 
of the analysis. In Section IIIII we discuss how possible 
complications arising from the use of P — > Kn decays 
as quasi-flavour-specific states can be resolved. In Sec- 
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tion IIVI wc describe the Dahtz plot models that we use 
in our feasibility study. We finally present the results of 
the study in Section |Vl before concluding in Section IVll 



II. FORMALISM 

We begin by recalling the essentials of the — > 
DK^, D -> A'^7r+7r~ model-independent method [l^, 
[il. The amplitude of the B+ ^ DK+, D if^Tr+Tr" 



decay can be written 



A 



_DDlz 



Ad + rBe*(*«+TUz 



where Ad = "i^^^-) = ADiml,m^_) is 

the amplitude of the — > A'^tt+tt" decay, Ad = 
^r>(™+,w^) is the amplitude of the D° — >• K^tt'^t:^ 
decay, is the ratio of the absolute values of the in- 
terfering i3+ — > D"K+ and i?+ — > D^K+ amplitudes, 
and Sb is the strong phase difference between these 
amplitudes. Assuming no CP violation in D decay, 
ADim\,'int) = AD{m'^,m\). The density of the D de- 
cay Dalitz plot from _B+ DK^ decay is given by the 
absolute value squared of the amplitude 



1^ 



DDlzl 



\AD+rBe'^'^+''^AD\^ = 



(2) 



+ rl\AD\^ + 2\AD\\AD\ixc - ys) , 



where 



X ^ rs cos{Sb + 7) ; y = rB sm{5B + 7) • (3) 
The functions c — ^^-y^'^ 0^,2 ^ „ _ ^f^2 



s(m^,7ni) are 



c(rn^ , iri^ ) and s 
the cosine and sine of the strong phase difference Sd ^ 
arg^D — arg^D between the — ATgTr+Tr" and Z?" — > 
KgTr~^TT~ amplitudes^: 

c = cos(5D(m^, m?_) ; s = sin 6 DiiTL^jin^) . (4) 

The equations for the charge-conjugate mode B^ — > 
DK~ are obtained with the substitution 7 — > —7. Con- 
sidering both B charges, one can obtain 7 and 5b sepa- 
rately. 

Once the Dalitz plot is divided into 2M bins symmet- 
rically to the exchange rrpL the expected number 
of events in the i**^ bin of the Dalitz plot of D decay from 
B+ DK+ is 



i_D Dlz 



Ki + r^K^t + 2^ K^K^i(xc^ - ysi) 



(5) 

where Ki is the number of events in the corresponding 
bin of the Dalitz plot where the D meson is in a flavour 



^ This paper follows the convention for strong phases in D decay 
amplitudes used by CLEO [20| . Note that alternative conven- 
tions are used in the literature, e.g. Ref. [iSj uses a convention 
where the sign of is opposite. 



eigenstate (conveniently obtained using D*^ -> Dtt^ 
samples) and /iddiz is a normalisation constant. The bin 
index i ranges from —M to M (excluding 0); the exchange 
rn?_ corresponds to the exchange i o —i. The co- 
efficients Ci and Si, which include information about the 
cosine and sine of the phase difference, are given by 



/ \Ad\\Ad \ cos 5 DdV 
J \AD\^dV ^ \Ad\HV 



(6) 



(1) M 



where the dependences on Dalitz plot position of \Ad\^ 
,Ad\ and 5d have been suppressed. The parameters 
Si are defined similarly with cosine substituted by sine. 
Here T) represents the Dalitz plot phase space and Vi is 
the bin region over which the integration is performed. 

The symmetry under tt"*" -f-)- 7r~ requires Ci = C-i and 
Si = — s_i, thus for 2M bins there are 2A/'+ 3 unknowns. 
For each bin there are two observables in B — DK de- 
cays, being the numbers of events in B^ and B~ samples, 
and thus there are a total of AJ\f observables. The system 
can consequently be solved for A/" > 2. However, due to 
the small value of tb in B^ — > DK^ decays, there is very 
little sensitivity to the Ci and Si parameters, leading to 
a reduction in the precision on 7 that can be obtained. 
Use of external constraints on Ci and Sj, that can be pro- 
vided by charm-factory experiments [20[, leads to much 
improved sensitivity to 7. 

Turning now to S° DK+TT~ decays, we find 
that the simplest approach towards performing a model- 
independent analysis is not viable. Using the convention 
Dcp = {D^ ± Z)") /V2, where -I- and — correspond to 
CP-even and CP-odd final states, respectively, the am- 
plitude for P° DcpK^TT~ decay is given by 



AbHIz — (±^S 



(7) 



where Ab = AB{iTijj^,m1^.^) is the amplitude of the 
P" -> WK+tt^ decay and Ab = AB{m%^,mj^.^) 
is the amplitude of the P" — > D^K'^tt~ decay. In 
these expressions we have factored out the CP- violating 
phase 7 so that the amplitudes Ab and Ab are CP- 
conserving. Both the strong phase difference between 
P° WK+t:- and P° D^K+n- decays and the ra- 
tio of the suppressed and allowed amplitudes vb are now 
functions of P decay Dalitz plot variables and enter the 
expressions for amplitudes. 

Let us now consider dividing the P decay Dalitz plot 
into bins (denoted by the index a, 1 < a < Ai). The 
number of expected events in each P'' — > DcpK^Tr~ 



Dalitz plot bin is 

(Ma) = l-lB Dlz 



2\J NaN a{xa COS 7 — (Tq, sin7) 



(8) 



In this expression the factors and N^, give the 
numbers of events in the corresponding bin in P" — >■ 
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D^K~^7T~ and — D'^K^tt~ decays, respectively 
(complications arising due to the difficulty to reconstruct 
D mesons in pure flavour eigenstates are discussed in 
Section IIIII below) . The factor h b diz is a global nor- 
malisation constant. In Eq. (|8]) we have introduced 
the parameters >Ca and CTq that describe the average 
strong phase difference and amplitude difference between 
B° WK+tt- and B° D^K+tt' decays in each bin, 
and that are defined by: 

/ \ Ab\\Ab\ cos SedV 

= 7° — = ' 

where the dependence on the B decay Dalitz plot posi- 
tion of \Ab\, \Ab\ and Sb = argA^ — arg^s have been 
suppressed, V again represents the phase space (now of 
the B decay Dalitz plot) and is the bin region over 
which the integration is performed. The parameters CTq, 
are given by similar expressions with cosine substituted 
by sine. 

Until this stage the discussion for i?" DK^tt^ 
decays has run in parallel to that for i?+ DK^, 
D Kl-K+TT^. However, unlike the B+ DK+, 
D — > ifgTT+TT" case, the two amplitudes Ab and Ab are 
inherently different and thus there is no symmetry in the 
Dalitz plot to be exploited. In the case of CP-eigenstate 
D decays, we therefore have 2Ai+l unknowns (even if the 
normalisation is fixed externally), that cannot be solved 
for with 2A4 observables. Hence external constraints on 
Ka and CTq, are necessary (l2| . 

Let us now consider both Dalitz plots simultaneously 
in the B° DK+tt~ , D -> A'^tt+tt" decay chain. The 
amplitude can be written as 

^DblDlz = 'Ab'Ad + c'^AbAd ■ (10) 

The decay probability is proportional to the absolute 
value squared of the amplitude A: 

lAoblDlzP - \Ab?\AdI + \Ab?\Ad? + 
2\Ab\\Ad\\Ab\\Ad\ X 

Re ^e*'^ e'*^ ''"o" ) ^^^d (m^ ) ^ 

= \Ab\'_\AdI + \Ab?\Ad? + 
2\Ab\\Ad\\Ab\\Ad\ X 

[{xc — as) cos 7 — {xs + ac) sin 7] . 

(11) 

To implement the model-independent analysis both B 
and D plots have to be binned. Corresponding bin in- 
dices are denoted here, as before, by Greek and Latin 
characters, respectively. The expected number of events 
in the bins of 5° DK+tt-, D KItt+tt- double 
Dalitz plot is then 

{Mai) = /lDblDlz{lV„A'j + NaK-, + 

2^NaIuNo,K^, X 

[{XaCi - aaSi) COS 7 - {XaSi + ffaCi) siu 7] | , 

(12) 



where the definition of the event numbers Na, Na, Ki 
and phase terms Ci^Si, Xa, a a is the same as for the pro- 
cesses mentioned above, and /lobiDiz is a normalisation 
constant. 

If the number of bins in the B^ DK^tt^ Dalitz 
plot is Ai and the number of bins in the D Kgir^n^ 
Dalitz plot is 2J\f, then after exploiting the symmetry of 
the D — )- KgTT'^Tr~ decay as before, the number of equa- 
tions represented by Eq. is M.M . The number of 
unknowns (including the normalisation factor ft-DbiDiz) is 
2M -\- 2Af + 2. Naturally, information on 7 cannot be 
extracted from the decays of only one B flavour, as is 
apparent since Eq. (|12p is invariant under the rotation of 
7 with the simultaneous rotation of all {x^ a) pairs. To 
determine the CP violation parameter it is, of course, 
necessary to add the opposite B flavour (with the substi- 
tution 7 — > —7). The system can now be fully resolved: 
in total, there are 2M.M observables and 2M. 4- 2M + 2 
unknowns, i.e. it is solvable for (A/ — 1)(A^ — 1) > 2. 
Since the coefficients ci and si are the same as in the 
model-independent B — > DK analysis, and can be ob- 
tained at a charm factory, the number of unknowns can 
be reduced even further. 

As a by-product of the i3° DK+t:^ , D A'^tt+tt" 
analysis, the coefficients x^ and ctq will be determined. 
The values of these parameters can then be used to en- 
able a model-independent analysis of B° — > DcpK^tt~ 
decays. Technically the combination of B'-^ DK^tt^, 
D — KgTr~^7T~ and — > DcpK^n^ modes is con- 
veniently done by using a combined likelihood fit that 
includes both the expressions Eq. (|12p and Eq. 



III. EFFECTS OF 
DOUBLY-CABIBBO-SUPPRESSED 
CONTRIBUTIONS TO D ^ Ktv DECAYS 

In practise it is hard to reconstruct secondary D 
mesons from B decays in pure flavour eigenstates. (The 
possibility to exploit the charge of the associated pion in 
£)*± _^ Dtt^ is not available.) Semileptonic D° decays in 
principle provide a source of pure flavour-specific states, 
but owing to small branching fractions and large back- 
grounds these are experimentally difficult to deal with, 
especially at hadron experiments such as LHCb. It is 
more convenient to use hadronic final states such as ATtt, 
which are not pure flavour tags due to the presence of 
a small admixture of doubly-Cabibbo-suppressed ampli- 
tude. This contribution can significantly affect a preci- 
sion measurement. 

In the proposed analysis, hadronic final states can be 
used after some correction of the procedure. The ampli- 
tude of a i?° DK+Tr~, D AT+tt" decay is 

^fav = Ab+ TKne-'^"- e'^ Ab , (13) 

where Vk-k and 6^^^ are the ratio of the magnitudes of 
the suppressed and favoured D decay amplitudes and the 
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strong phase between them, respectively.^ The number 
of events in the B'^ — > DK^tt^ Dahtz plot bins with D 
detected in a K^n~ state is: 



coa{-6KTv + l) - (^a sin(-(5if7r + 7)] ■ 

(14) 

Similarly for the suppressed decay B'^ — > DK^tt , 
D — > K^TT^, the amplitude is 



and the number of events is 



(15) 



[>fa COs((5K7r + 7) - sin((5A'7r + 7)] ■ 

(16) 

The expressions for B^ DK 7r+ decays are obtained 
by the substitution 7 — > —7. 

The relations of Eq. ([H]) and Eq. add another AM 
equations to the system of equations to be used in the 
analysis {2M equations for each B flavour), but since 
the constants tk-k and 5kit are known from charm anal- 
yses [111, [l^l, only 27W new unknowns are added (the 
numbers and No). Hence, the inclusion of the sup- 
pressed final states into the analysis can help to improve 
the overall sensitivity to 7 [2^ [l^] • Recent studies sug- 
gest further improvement could be gained from including 
also the D — > Ktttt'^ decays 25 1. 



IV. B DECAY DALITZ PLOT MODELS 

In order to perform a feasibility study of the proposed 
method it is necessary to define a B decay Dalitz plot 
model. The sensitivity depends crucially on the struc- 
ture of the B° D^K+TT- and B° D"K+tt- ampli- 
tudes, specifically on the interference between them. If 
the strong phase difference Sb is nearly constant {e.g. 
in the case when both amplitudes are dominated by 
B — > DK* decay), the double Dalitz analysis reduces 
to a conventional Dalitz plot analysis technique similar 
that used for B^ DK^ decays. Note that in this 
type of analysis the possible contribution from non-A'* 
states can be taken into account model-independently by 
introducing the coherence factor < « < 1 to the inter- 
ference term [l3| (analyses using this approach have been 
performed, albeit with low statistics 0, [131). However, 
if the resonance structure in B° — )- DK'^-k~ decays is 
rich with multiple overlapping states in different chan- 
nels [K* , D* and £**), the proposed technique should 
provide a significant benefit compared to analysis of the 
D — >■ iCoTr+Tr" Dalitz plot with the coherence factor. 



^ We use a sign convention for Skti that is consistent with that 
used in the majority of the literature (see, e.g. Ref. but is 
opposite in sign to that used for Su . Another different convention 
is used by CLEO [HI,!!!. 



Unlike the D — A'^Tr+Tr" amplitude, the amplitudes 
for B^ — )• DK~^TT~ decays are not yet well established. 
For our feasibility study, based on pseudo-experiments 
generated using Monte Carlo (MC) simulation, we use 
several models based on two-body amplitudes with K* 
and D* resonances. We define first model A, which in- 
cludes only A'* (892) and £'2(2460) states in the favoured 
amplitude and K*{892) in the suppressed amplitude with 
relative magnitudes consistent with measured values for 
the branching fractions [2^ and phases chosen arbitrar- 
ily. Model B is similar to that used in a recent study [O] , 
where the contributions of further resonances are based 
on results from B — >■ D^K^ Tr~ [2^ using the assump- 
tion of isotopic symmetry. Model C includes a contribu- 
tion from the Ds.j (2700) resonance [s^l in addition to all 
states of model B. Its contribution is estimated from the 
ratio 



B{B Ds.jn) _ B{B Di*\) 
B{B^D,jD) ^ B{B^Di*^D) 



(17) 



with a phase space correction. Models D and E use the 
same states as models B and C, respectively, excluding 
the nonresonant term. For the D — >■ KgTr~^TT~ amplitude, 
we use the model from the most recent Belle analysis [U . 

The fit fractions and phases for model C are shown 
in Table HI Other models contain subsets of the states 
from model C with the same amplitudes and phases. The 
Dalitz plots of 5° ^ WK+tt' and 5° ^ D^K+n' de- 
cays generated according to model C are shown in Fig- 
ure [TJ 

The sensitivity of the technique will depend on the 
choice of binning for both D — >■ Kgir^ir^ and i?" 
DK~^7r~ Dalitz plots. The choice of optimal binning for 
D — > KgTT^TT^ Dalitz plot has been previously consid- 
ered in application to B~^ — >■ DK~^ modes [13. Here 
the requirements for the optimal binning are the same: 
the strong phase difference {6b in case of — >■ DK~^tt~ 
plot and 6d for D — >■ Kg7T~^TT~ plot) and the ratio of 
suppressed to allowed amplitudes should be as constant 
as possible over each bin to maximise values of x'^ + cr^ 
and cf + sf. In the case of the D — KgTr~^7T~ Dalitz plot, 
a binning based on uniform division of strong phase dif- 
ference has been found to give a good approximation to 
the optimal [l^. However, in the case of B^ — >■ DK~^Tr~ 
amplitudes, such an approach gives poor results because 
of very large variations of the absolute value of amplitude 
ratio across the bin. A binning based on maximising the 
"binning quality factor" Q [l3|, gives significantly bet- 
ter results. The optimal binning of the i?" DK'^7t~ 
Dalitz plot with 7 bins for model C is shown in Figure [2] 
(a), while Figure [2] (b) shows the values of {xa,(Ja) ob- 
tained in each bin. For the D ATgTr+Tr" Dalitz plot, we 
also use optimised binning obtained with the same tech- 
nique (see Figure [3]) with binning quality factor Q=0.89. 

We use the values vb ~ 0.4 and 7 = 60° in the sim- 
ulation. Note that the value of can differ for differ- 
ent K* states {e.g. K*{892), iq{U30), A:2*(1430) etc.); 
we use the same value for all these states. We generate 
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FIG. 1: Generated Dalitz plots for (a) B° D^A'+tt" and (b) B° D^K+tv' decays using amplitude model C. 
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FIG. 2: (a) Binning of — ^ DK^tt Dalitz plot and (b) values of Xi, ai parameters for optimal binning (crosses are calculated 
values, scattered points are the values obtained from the fit to toy MC samples). 




m^K.JGeV^/c") <=: 

FIG. 3: (a) Binning of D ^ K%-k'^'k~ Dalitz plot and (b) values of d, Si parameters for optimal binning (crosses are calculated 
values, scattered points are the values obtained from the fit to toy MC samples). 
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TABLE I: Parameters of B° — > DK^% amplitudes. The suppressed B° D^K^-k amplitude uses rs — 0.4 





B" ^ 


D"A'+ 


TT 








Amplitude 


Phase 


Fit frac. 


Amplitude 


Phase 


Fit frac. 


7<'*(892) 


1 


0° 


47% 


rB X 1 


120° 


56% 


i^o*(1430) 


0.2 


284° 


0.015% 


TB X 0.2 


284° 


0.018% 


7^2(1430) 


0.16 


221° 


10% 


tb X 0.16 


221° 


12% 


7^*(1680) 


0.61 


128° 


2.4% 


rs X 0.61 


128° 


2.9% 


L»S(2400) 


4.8 


267° 


2.3% 








L»2(2460) 


1.1 


325° 


30% 








Dsj(2700) 








TB X 4.0 


200° 


25% 


Non-res 


1.25 


140° 


5.3% 


r-B X 1.25 


140° 


6.3% 



pseudo-experiments with samples of 10^ events for each B 
flavour for both DK+n^ with D /irgTr+Tr" and 

— )■ DcpK^-K" modes. Note that once the Dahtz plot 
models are fixed, the ratios of the numbers of events be- 
tween the different samples is also fixed. The numbers of 
B° WK+tt- and B^ D^K+tt' decays are 10^ and 
~ 8000, respectively. These numbers of events roughly 
correspond to 20 times the expected annual yield of the 
LHCb experiment.^ Hence the precision obtained from 
this study is indicative of that which could be achieved 
in an upgraded phase of LHCb operation. It is, however, 
trivial to scale the sensitivities to lower luminosities. 



V. RESULTS OF THE FEASIBILITY STUDY 

Different strategies for using B'^ — > DK^n^ decays to 
measure 7 are studied using MC simulation. The results 
are shown in Table [Hi The uncertainties are obtained 
from the spread of results in the samples of 200 pseudo- 
experiments. 

An unbinned fit of — > DcpK^n^ decays with free 
parameters r^, 7 and 5b (option 1) gives an estimate 
of the statistical error for a model-dependent Dalitz plot 
analysis of this decay . A binned fit of the same mode 
(option 2) with fixed values of and a shows a drop of 
sensitivity due to binning of the — > DK^tt^ Dalitz 
plot. For most amplitude models, the loss in sensitivity 
is about 20%. 

We next consider selecting only events in the K* region 
of i?° — > DK^TT~ phase space, introducing a coherence 
factor [l^l, and using D — > A'^tt+tt^ decays. A conven- 
tional single Dalitz plot analysis (option 3) gives com- 
parable sensitivity to the unbinned B^ — >• DcpK^n^ 



* A full simulation of the detector response to — )■ DK~^it~ , 
D — > A^TT+TT" decays has not been performed, but we assume 
that the ratio of the efficiency of this mode to — > DcpK^ir^ 
mode is equal to that between B"*" — > DK'^ , D — >■ Kg7T^n~ and 
i?"*" — > DnpK^ decays, where roughly equal annual yields are 
expected |32|| . 



analysis (option 1). However, only a fraction of all 
B^ — > DK^TT^ events are used indicating that further 
gains are potentially possible. Changing from unbinned 
to a binned model-independent approach (option 4) re- 
duces the sensitivity by about 15%. Relaxing the Ci and 
Si coefficients in the binned approach is also possible (op- 
tion 5) Although the sensitivity in this case re- 
duces significantly, the reduction is far less than in the 
B^ —J- DK^ mode [s^ due to the larger interference 
term in neutral B decays and also due to the optimal 
binning used in this study. 

The double Dalitz plot analysis (option 6) uses all 
available events in the i?" — > DK^n~ phase space, pro- 
viding better 7 accuracy compared to option 4, and also 
giving as a by-product the values of B^ — >■ DK'^Tr~ am- 
plitude coefficients >Ca and CTq,. The values of and ctq 
returned from the fit and a comparison with their cal- 
culated values is shown in Figure [5] (b) . Relaxing the 
D — > KgTT^TT^ amplitude coefficients q, s; (option 7) re- 
sults in a much smaller reduction of the sensitivity than 
that with single Dalitz plot analysis selecting only the 
K* region. The fitted values of the Ci, Si coefficients for 
this option arc shown in Figure [3] (b). 

Combined with B'^ — > DcpK^ir^ decays (option 8), 
the double Dalitz analysis allows to further improve the 
statistical precision. However the gain in sensitivity for 
most models is about 15% worse than the weighted error 
for options 6 and 2, which means that the precision of Xa 
and ctq, terms obtained from the Dalitz analysis affects 
the B^ — > DcpK^TT^ mode. With very large statistics 
it may be possible to mitigate this effect by using a finer 
binning of the B° DK^tt~ Dalitz plot. 

An example of the residual distribution for 7 using 
option 6 (double Dalitz plot fit) with amplitude model 
C is shown in Figure |3K a). No systematic bias is seen in 
any of the fits. 

Wc emphasise that, although the fit procedure uses 
a binning that depends on the decay amplitude, the 
method does not generate systematic bias even if the am- 
plitude used to define the binning is wrong or the model- 
blind binning is used. To demonstrate this, we perform 
the fit of the data generated using model C, but using 
model-blind rectangular Dalitz plot bins (the bin bound- 
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TABLE II: Sensitivity of 7 measurement with different techniques using B° decays 







Model A 


Model B 


Model C 


Model D 


Model E 




Composition 


K*, 
D2 only 


No Dsj 


All states 
included 


No Dsj and 
no nonres. 


No nonres. 




_B" DK^-K~ binning quality Q 


0.90 


0.91 


0.81 


0.86 


0.77 


1 


Unbinncd fit of B" DcpK+ix- 


2.47° ±0.16° 


1.31° ±0.07° 


1.17° ±0.06° 


1.71° ±0.09° 


1.53° ± 0.08° 


2 


Binned fit of B° DcpK+n- 
(fixed yc, (j) 


3.45° ± 0.33° 


1.56° ±0.08° 


1.46° ± 0.07° 


2.07° ±0.11° 


1.86° ± 0.10° 


3 


Unbinned Dalitz fit of B° DK*°, 

D K^TT+TT- 


1.34° ± 0.07° 


1.50° ± 0.08° 


1.54° ±0.08° 


1.46° ±0.08° 


1.50° ±0.08° 


4 


Binned Dalitz fit of B" DK'°, 

D K^TT + TV- 


1.48° ±0.07° 


1.71° ±0.09° 


1.73° ±0.09° 


1.72° ±0.09° 


1.78° ±0.09° 


5 


Binned B° DK'^, D J^gyr+Tr^ fit 
(floated Ci, Si) 


2.03° ± 0.11° 


2.67° ± 0.18° 


2.54° ±0.16° 


2.40° ±0.15° 


2.40° ±0.15° 


6 


Binned double Dalitz fit 
(fioated >c, a) 


1.37° ± 0.09° 


1.49° ± 0.08° 


1.45° ± 0.07° 


1.46° ±0.07° 


1.50° ± 0.08° 


7 


Binned double Dalitz fit 
(fioated >c, a, Ci, Si) 


1.81° ±0.09° 


1.56° ± 0.08° 


1.51° ±0.08° 


1.51° ±0.08° 


1.59° ±0.08° 


8 


Binned double Dalitz fit with 

DcpK+n- (floated K,a) 


1.29° ±0.06° 


1.21° ±0.06° 


1.20° ± 0.06° 


1.39° ±0.07° 


1.31° ± 0.07° 


9 


Binned double Dalitz fit with 

B" — >• DcpK^Tv^ (fioated d, Si) 


1.52° ±0.08° 


1.27° ±0.06° 


1.24° ± 0.06° 


1.36° ±0.07° 


1.41° ±0.07° 



aries in mj^^ arc ]\/P{K* (892)) and 7 McV^/c'^; those in 
mjj^ are AP{D^{2i60)) and 12 MeVVc^ the bins above 
"^Ktt = MeV^/c^ are combined so that the number 
of bins is 7 as in the other fits). The quality factor for 
this binning is Q = 0.60. The residual distribution for 7 
using option 6 (double Dalitz plot fit) is shown in Fig- 
ure |4l^b). The fit is unbiased, and the statistical accuracy 
compared to the fit with binning based on the correct un- 
derlying model is roughly proportional to the ratio of the 
Q values: the 7 resolution is (7(7) = 1.93 ± 0.10°. 

The statistical accuracy of the method does not depend 
significantly on the value of 7. The precision on 7 as a 
function of 7 is shown in Figure SI^c) for model C, options 
6 (double Dalitz plot fit only) and 8 (double Dalitz plot 
fit combined with i?° — > DcpK~^tt~). We expect that 
the dependence on the strong phases in B'^ DK^tt^ 
decays, which were taken arbitrarily in our study, should 
also not be large. 

VI. CONCLUSION 

We have proposed a technique to obtain the angle 7 
of the Unitarity Triangle in a model-independent way by 
using the decay chain — >■ DK~^Tr~, D — >■ KgTT~^Tr~. 
Our method is independent of model assumptions in 



both B° — > DA'+TT^ and D — > /y^tt+tt^ amplitudes. 
The proposed approach allows not only to measure 7 in 
this decay chain, but also to use the constraints on the 
B'^ — > DK~^'!T~ amplitude obtained in this measurement 
for a model-independent 5° DcpK^n^ Dalitz plot 
analysis. 

The precision of the 7 measurement using our tech- 
nique will depend strongly on the structure of 5° 
DK^TT~ amplitude, which is as-yet unknown. However, 
MC simulation results show that if interference effects 
between B° D^K+tt' and B° WK+n' ampli- 
tudes are strong enough to allow a meaningful model- 
dependent measurement of 7 using B'^ — > DcpK^tt^ de- 
cays , and the number of reconstructed B° — > DK^n^ , 
D KgTr'^7T~ decays is comparable to the number of 
B'^ — > Dc'pK'^Tr~ decays, the model independent tech- 
nique proposed here should enable a significant improve- 
ment in the sensitivity. 
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and Mike Williams for discussions. We would like to ac- 
knowledge support from the University of Warwick Re- 
search Development Fund, the Science and Technology 
Facilities Council (United Kingdom) and from the Min- 
istry of Science and Education of the Russian Federation. 
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